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Abstract: The Greenberger-Horne-Zeilinger (GHZ) effect provides an example of quantum corre- 
lations that cannot be explained by classical local hidden variables. This paper reports on the 
experimental realization of GHZ correlations using nuclear magnetic resonance (NMR). The NMR 
experiment differs from the originally proposed GHZ experiment in several ways: it is performed 

. on mixed states rather than pure states; and instead of being widely separated, the spins on which 

it is performed are all located in the same molecule. As a result, the NMR version of the GHZ 
experiment cannot entirely rule out classical local hidden variables. It nonetheless provides an 
unambiguous demonstration of the "paradoxical" GHZ correlations, and shows that any classi- 
cal hidden variables must communicate by non-standard and previously undetected forces. The 
NMR demonstration of GHZ correlations shows the power of NMR quantum information processing 

■ techniques for demonstrating fundamental effects in quantum mechanics. 

Quantum mechanics is well-known to exhibit strange and apparently paradoxical effects. Per- 
haps the strangest of these effects is the possibility of correlations between quantum-mechanical 
systems that cannot be explained by appealing to an underlying theory of local hidden variables — 
classical variables whose unknown or "hidden" states supply the apparent stochastic results to 
quantum experiments ^ |], ||, A striking contradiction with the local hidden variable 

J> | picture is obtained in the Greenberger-Horne-Zeilinger (GHZ) effect, in which a set of measure- 

ments on three correlated quantum particles gives results that are perfectly correlated in a pattern 
that is classically impossible |7|. This paper reports the first experimental demonstration of GHZ 
correlations. 

GHZ experiments that verify both the "weird" quantum correlations and the non-local nature of 
physical reality are difficult to realize in the laboratory, where the problems associated with faith- 
fully preparing the required state \ip) and making the necessary measurements pose formidable 
obstacles. However, recent discoveries in the field of quantum information processing using NMR 
make possible an analog GHZ experiment that validates the GHZ correlations required by a quan- 
tum mechanical description of the universe H,|, 0. The idea of using NMR to demonstrate GHZ 
correlations was proposed in ]ll|, and an effective GHZ state was first created using NMR in fl2f . 
NMR offers a number of advantages in performing tests of few-particle quantum mechanics, in- 
cluding: (1) weak measurements that allow the direct determination of ensemble density matrices; 
(2) long decoherence times; (3) simple, well-understood Hamiltonians; (4) gradient techniques that 
allow the construction of effective pure states and the controlled decoherence of samples. These 
features of NMR make it possible to use quantum systems in mixed states at room temperature to 
mimic precisely the dynamics of a quantum system in a pure state at low temperature. As a result, 
NMR is an ideal testbed for the predictions of quantum mechanics for few-variable systems. 
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Despite the advantages of NMR quantum information processing techniques listed above, it is 
important to note that demonstrating GHZ correlations using NMR differs in significant ways from 
performing the GHZ experiment as originally envisaged. In particular, because the NMR experi- 
ment is inherently local, the NMR demonstration of GHZ correlations cannot rule out nonstandard 
interactions between classical hidden variables that conspire to exhibit apparently quantum me- 
chanical correlations. In addition, because the NMR experiment is performed on thermal states, 
though the spins in the experiment mimic the effects of entanglement, they are not in fact entan- 
gled: as a result, the NMR experiment cannot rule out local hidden variables. Nonetheless, the 
NMR demonstration of GHZ correlations provides an unambiguous and statistically significant test 
of the predictions of quantum mechanics. The experiment shows these to be accurate to a high 
degree of precision. 



The following version of the GHZ experiment is due to Mermin [13]. Consider three spin- 1/2 
particles prepared in the state 

1^) = ^ (]+ + +) - | » , (1) 

where |+) represents "up" along the z-axis and |— ) signifies "down" along the z-axis. Now consider 
the results of the following products of spin measurements, each made on state 

(i) particle (1) along x, particle (2) along y, and particle (3) along y. Note that since a\a^ j o^\)) = 
+1^)) the product of the results of these measurements should be +1, where +1 corresponds to the 
result |+), and —1 corresponds to |— ). 

(ii) particle (1) along y, particle (2) along x, and particle (3) along y. As in scenario (i), since 
a l a x cr y\ l l') = +1^)] the product of these results should be +1. 

(iii) particle (1) along y, particle (2) along y, and particle (3) along x. Again, since cr^o^o"^ [-(/>) = 
+\tp), the product of these results should also be +1. 

(iv) particle (1) along x, particle (2) along x, and particle (3) along x. In this case, there is a 
significant sign difference: a^.a^.a^l'ip} = so the product of the results should be —1. 

The minus sign in the final scenario is crucial for differentiating between quantum mechanical 
and classical (i.e. hidden- variable) descriptions of reality. No classical hidden- variable model that 
assigns a classical variable (such as +1 or —1) to the measurement of each individual particle's spin 
along axes x and y could reproduce the quantum predictions: classically, the product of the four 
scenarios above for a hidden- variable model must yield +1 because each particle is measured along 
each of the two axes exactly twice. 

In addition to demonstrating purely quantum correlations that have no classical analog, the 
above set of measurements has the potential to verify the hallmark of quantum "weirdness" : non- 
locality. The GHZ experiment — as originally conceived — requires that the measuring apparatuses 
for the three particles be sufficiently distant from each other so as to eliminate the possibility of 
local, classical interactions between measurements that might duplicate the quantum mechanical 
predictions. In NMR experiments, the effective distances between any three spins in a molecule 
are typically a few angstroms at most; as such, it is practically impossible to demonstrate non- 
local effects using NMR. However, in the NMR experiment below, the system's Hamiltonian is well 
known to a high degree of accuracy; thus, the NMR experiment can rule out any standard local 
interaction as the cause of the non-classical correlations. 

This paper reports on the first experimental demonstration of GHZ correlations using effective 
pure states. Techniques now exist whereby it is possible to create an effective pure state in a liquid 
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sample of 0(Na) essentially non-interacting molecules at room temperature, meaning that the 
state of the ensemble can be described by a density matrix whose deviation from unity transforms 
under unitary operations exactly as a pure state density matrix. Since only the deviation density 
matrix (p), averaged over the whole sample, contains observable magnetization, quantum mechanics 
predicts that the results obtained from performing unitary transformations on a pseudo-pure NMR 
state are identical to those that would be found were one to perform the same operations on a single 
quantum system. For example, the same unitary operator (e.g., IA = e % ^ ax(Txa y) that transforms a 
pure state (e.g., |+ + +)) into a GHZ state (e.g., also transforms the deviation density matrix 
of an ensemble pseudo-pure state into a density matrix with the identical, non-classical correlations. 

It is important to keep in mind that although the density matrix for an ensemble in an effective 
pure GHZ state is the same as the density matrix for an ensemble in which some fraction of the 
systems are in a pure GHZ state and the remainder are in a completely mixed state, these two 
ensembles are not the same. Two different ensembles may have the same density matrix. In the 
ensemble constructed here, each of the systems is in a unitarily-transformed version of a thermal 
state; none are in GHZ states. Nevertheless, since two ensembles with the same density matrix 
are indistinguishable with respect to measurement, the effective GHZ state constructed is just as 
effective as an actual GHZ state for testing the predictions of quantum mechanics. 

Our verification of GHZ correlations was performed on a Bruker AMX400 spectrometer at room 
temperature. At this temperature, the Boltzmann distribution over an ensemble of homonuclear 
spins in a liquid NMR sample gives the average density matrix (per molecule) approximately 
proportional to 1 + J2j fij^i, where (3j = 2pjB/ksT = hujj/kBT rs huj/ksT w (3 is the Boltzmann 
factor for the j'th spin. With uj ~ 100 MHz for 13 C in a static magnetic field of B = 9.6 Tesla, the 
room-temperature (T ~ 300 K) Boltzmann factor is of order 10~ 6 . Higher order terms in the density 
matrix expansion are at least of order 10~ 12 , and can therefore be neglected. Only the deviation 
of the density matrix from unity, rescaled in the rest of this paper (at equilibrium, p eq = J2 a i)i 
represents surplus or deficit populations in energy levels whose transitions are observable in NMR. 

In our experimental setup, one set of coils was available to provide a gradient in B z across the 
sample. The result of a gradient is to introduce phase variations in the sample that vary spatially, 
which, when averaged over the whole sample, represent non-unitary operations on the ensemble 
density matrix. In addition, another set of coils in the x — y plane can provide rf pulses that cause 
nearly unitary rotations of the density matrix. The rf coils are also used as pick-up coils to measure 
magnetization in the x — y plane. For three spins (as in the demonstration below), the observable 
magnetization is the spatially-averaged signal during NMR signal acquisition that is proportional to 



spins can be read directly from the spectra produced by plotting the Fourier transform of the 
induction signal. 

The sample used in the demonstration of GHZ correlations was triply-labled ( 13 C) alanine with 
X H decoupling. The lowest-order terms in the natural Hamiltonian for the system are 



where ujj is the Larmor precession frequency (including the effects of chemical shift) for spin j, and 



interactions between the carbon spins in alanine are Jyi = 53.4 Hz and J23 = 35.3 Hz. The scalar 
interaction between spins 1 and 3 is small (J 13 = 1.4 Hz); this term and all other interactions in the 
Hamiltonian (including the effects of diffusion and relaxation) can be neglected over the « 160 ms 
during which the entire experiment takes place. 




Other elements of the density matrix are unobserved. Correlations between 
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The experiment works as follows: 



(A) Equilibration. The sample is allowed to come to thermal equilibrium, in which the density 
matrix is: 



Peq 



^! _l ^ 2 , 3 

O?. + ^7. + °7 



( 3 




















\ 





1 


























1 


























-1 


























1 


























-1 


























-1 





\o 




















-3/ 



(3) 



(B) Preparation. The sample is prepared in an incoherent mixture of two pseudopure states 
using the following operations: 
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In this notation, each bracketed expression represents one matrix operation. Brackets containing 
angles indicate an rf pulse that selectively rotates the superscripted spin by the specified angle, 
where the subscripted "phase" designates the axis of rotation in the x — y plane of the co-rotating 
frame. Selective pulses in this experiment required a pulse length of 2 ms. Brackets with partial 
derivatives indicate an applied gradient in the magnetic field along the direction in the denominator. 
Brackets containing a coupling constant indicate free evolution of the system under the natural 
Hamiltonian for the time indicated; this evolution allows the pertinent scalar coupling term to 
correlate spins. Pulses for refocussing or for averaging to zero other coupling terms in the natural 
Hamiltonian are not shown. 

The result of the preparation sequence is the density matrix: 
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which transforms like a balanced mixture of two pseudopure states, 
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(C) Rotation. The density matrix for the sample is now rotated via unitary operator e 1 a^^^. 
using the sequence: 
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The density matrix thus rotated becomes 
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Note that this matrix has the identical off-diagonal structure as the density matrix formed from 
the pure state of Eq. |l|: 

1 



Pghz 
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i.e., the NMR sample as prepared above represents the same quantum correlations that would be 
present in a pure, three-particle GHZ state. 

(D) Measurement. An NMR measurement that demonstrates one of the {a\cry(T^\ , {ayO^a^j , 

(ayayU^) , or (o-^a^a^.) correlations is now performed by applying rf pulses to the sample to rotate 

the desired correlation into observable magnetization. Specifically, to measure fcrj<j|of ^ , the pulse 
sequence 
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followed by data acquisition with the absorptive signal phased along k will reveal the magnitude and 
sign of the desired correlation in the spectrum on resonance with the second spin. Steps (A)~(C) 
of the experiment are then identically repeated three times, each with a different measurement 
sequence (D). 

The results of the four measurements on the sample prepared in state pghz are shown in Fig. 1. 
For example, to determine the correlations between {a^o^a^j, look at Fig. 1(a): the spectrum 
shows a multiplet of four lines. The frequency of each of the four lines corresponds to the resonant 
frequency of the second spin when the first and third spins are, respectively (from left to right): 
I - ) 3 > |+)i| — ) 3 > | — )i|+) 3) and |+)i|+) 3 . The fact that the phase of the first line from the left 
is (+) confirms the GHZ correlation that when spin 1 is |— ) along x and spin 3 is |— ) along y, spin 

2 is |+) along y. Similarly, the second line (from the left) shows a negative peak, confirming the 
correlation that when spin 1 is |+) along x and spin 3 is |— ) along y, spin 2 is |— ) along y. Likewise, 
the third line, with a negative peak, confirms the GHZ correlation that when spin 1 is |— ) along x 
and spin 3 is |+) along y, spin 2 is |— ) along y. Finally, the last line confirms the correlation that 
when spin 1 is |+) along x and spin 3 is |+) along y, spin 2 is |+) along y. Thus all four lines in 
Fig. 1(a) confirm that the product of the spin of particle 1 along x, particle 2 along y, and particle 

3 along y yields the result +1, as in scenario (i) of the GHZ experiment explained above. 
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Demonstration of GHZ Correlations 
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Figure 1: The confirmation of the quantum prediction of GHZ correlations. The four figures show 
the spectrum of spin 2 for each of the measurements. Each spectrum is split into a multiplet of 
four lines, corresponding to the resonant frequencies for spin 2 when spins 1 and 3 are in the states 
| — )3, |+)i|— )3j | — )i|+)3, and |+)i|+)3 respectively. An 'up' line indicates the state |+) for 
spin 2, and a 'down' line indicates the state |— ). Plot (a) confirms the GHZ correlation that the 
product of measurements of a x on the first spin, a y on the second spin, and a y on the third spin 
yields the result +1. For example, the first line in (a) indicates that when spins 1 and 3 give the 
result (— l)i(— 1)3, spin 2 gives the result +1; the second line indicates that when spins 1 and 3 give 
the result (+l)i(— 1)3, spin 2 gives the result —1; etc., so that the product of the results is always 
+1. Similarly, (b) and (c) confirm that the products a^a^ay and a^a^a^ also yield the result +1. 
Plot (d), by contrast, shows that the product cr^a^a^. yields the result —1, in contradiction to the 
predictions of the classical hidden variable theory. 
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Note that Fig. 1(b) shows the same four lines in the spectrum for spin 2, except that in this case 
the measurement sequence reveals the correlation that corresponds to a^a^ay. As in Fig. 1(a), all 
four lines indicate that the product of the spin of particle 1 along y, particle 2 along x, and particle 
3 along y is +1, as predicted in scenario (ii) above. Likewise, Fig. 1(c) shows the analogous +1 
result when the product of spin measurements of particle 1 along y, particle 2 along y, and particle 
3 along x is computed. 

Significantly, Fig. 1(d) shows a different spectrum when a measurement of the correlation 
{o\.o 2 x o 2 ^) was made. In this case, when spin 1 is |— ) along x and spin 3 is |— ) along x, spin 
2 is |— ) along x, as seen in the first line; when spin 1 is |+) along x and spin 3 is |— ) along x, 
spin 2 is |+) along x; and so forth. In this case, the product of the spin measurements of particle 
1 along x, particle 2 along x, and particle 3 along x is not +1, but rather —1. Taken together, the 
four spectra in Fig. 1 clearly demonstrate the GHZ correlations predicted by quantum mechanics. 

Although the NMR experiment demonstrates GHZ correlations, it cannot eliminate local in- 
teractions as the cause of those correlations because of the close proximity of the nuclei involved. 
However, no known interaction in the natural Hamiltonian could be responsible for the "communi- 
cation" between the nuclei that caused the GHZ correlations. Specifically, measurement step (D) 
required only 4.043 ms to implement; this represents the amount of time after the GHZ prepara- 
tion steps (A)~(C) — common to all four experiments — but before data acquisition, in which the 
culpable classical local interaction could "inform" the spins of which axes were being measured for 
which spins and then adjust the phase of the acquisition signal to reveal the quantum mechanical 
correlations. Since the fastest coupling term in the natural Hamiltonian is slower than 4.043 ms 
(I/2J12 = 9.36ms), the culpable "classical demon," if it exists, would have to be a non-standard 
interaction. 

In addition, the NMR experiment cannot completely rule out classical hidden variables. The 
effective pure state used to demonstrate GHZ correlations is in fact a mixed state of an ensemble 
in which individual members are unitarily transformed from a thermal state. It can be shown that 
such states can be described by non-negative discrete Wigner functions [16, 17]. Accordingly, the 
results of the experiment could be "explained" by a hypothetical ensemble of classical systems each 
of which possesses a definite state for its hidden variables. As in the case of local interactions in 
the previous paragraph, however, to be consistent with the measured properties of the states (3-5) 
above the preparation of such an ensemble from the thermal ensemble with which the experiment 
begins would require a conspiracy of non-standard interactions. 

In conclusion, this paper reported on an experiment that displays GHZ correlations using nu- 
clear magnetic resonance. Although the NMR experiment performed cannot entirely rule out 
classical local hidden variables, it nonetheless provides an unambiguous and statistically signifi- 
cant demonstration of the apparently paradoxical GHZ correlations. In addition, the experimental 
demonstration of GHZ correlations using NMR underscores the ability of NMR to perform signifi- 
cant experiments in fundamental quantum mechanics. 
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